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ABSTRACT 

We give an explanation for the origin of various properties observed in local infrared galaxies, and make pre- 
dictions for galaxy counts and cosmic background radiation (CBR), by a new model extended from that for 
optical/near-infrared galaxies. Important new characteristics of this study are that (1) mass scale dependence of 
dust extinction is introduced based on the size-luminosity relation of optical galaxies, and that (2) the big grain 
dust temperature rdust is calculated based on a physical consideration for energy balance, rather than using the em- 
pirical relation between Tdust and total infrared luminosity Lir found in local galaxies, which has been employed 
in most of previous works. Consequently, the local properties of infrared galaxies, i.e., optical/infrared luminosity 
ratios, LiR-Tdust correlation, and infrared luminosity function are outputs predicted by the model, while these have 
been inputs in a number of previous models. Our model indeed reproduces these local properties reasonably well. 
Then we make predictions for faint infrared counts (in 15, 60, 90, 170, 450, and 850 /im) and CBR by this model. 
We found considerably different results from most of previous works based on the empirical LiR-7dust relation; 
especially, it is shown that the dust temperature of starbursting primordial elliptical galaxies is expected to be 
very high (40-80K), as often seen in starburst galaxies or ultra luminous infrared galaxies in the local and high-z 
universe. This indicates that intense starbursts of forming elliptical galaxies should have occurred at z ^ 2-3, 
in contrast to the previous results that significant starbursts beyond z ^ I tend to overproduce the far-infrared 
(FIR) CBR detected by COBE/FIRAS. On the other hand, our model predicts that the mid-infrared (MIR) flux 
from warm/nonequilibrium dust is relatively weak in such galaxies making FIR CBR, and this effect reconciles 
the prima facie conflict between the upper limit on MIR CBR from TeV gamma-ray observations and the COBE 
detections of FIR CBR. The intergalactic optical depth of TeV gamma-rays based on our model is also presented. 

Subject headings: cosmology: observations — galaxies: evolution — galaxies: formation 



1. INTRODUCTION 

To understand when and how stars and galaxies formed in 
the universe is one of the most fundamental issues in modern 
astronomy. The energy emitted from stars as a result of nuclear 
fusion is radiated in two modes: direct stellar emission rang- 

. ing from optical to near-infrared (NIR) wavelengths, and emis- 
sion in mid- and far-infrared (MIR and FIR) wavelengths from 

, dust particles heated by stellar radiation field. ^ Faint galaxy 

■ counts and the cosmic background radiation (CBR), as well as 

■ the properties of galaxies in the local universe, give us impor- 
tant clues to understand galaxy formation. 

In optical/NIR wavelengths, the sensitivity of existing tele- 
scopes now reaches a depth sufficient to resolve more than 80- 
90% of CBR from galaxies (Totani et al. 2001a), thanks to very 
deep optical surveys such as the Hubble Deep Fields (HDFs: 
Williams et al. 1996; Williams et al. 2000) and the NIR sur- 
veys such as the Subaru Deep Field (SDF: Maihara et al. 2001). 
There are a number of theoretical models in various approaches 
to be compared with these data, and in fact these data can be 
reasonably explained in the framework of the Big-Bang cos- 
mology and structure formation induced by the cold dark matter 



(CDM). In contrast, infrared observations have not yet reached 
such depth and star formation activity hidden by dust extinction 
is still poorly known, despite the dramatic progress of obser- 
vations achieved by satellites such as IRAS (e.g., Soifer et al. 
1987 and references therein) and ISO (e.g. Puget et al. 1999; 
Oliver et al. 2000; Okuda 2000; for a review, see e.g., Gen- 
zel & Cesarsky 2000). The SCUBA of JCMT (Holland et al. 
1999) has opened a new window of submillimeter wavelengths 
to probe the dusty galaxies at high-z, and it seems that a con- 
siderable part of CBR from submm galaxies has been resolved 
(e.g.. Small, Ivison, & Blain 1997; Hughes et al. 1998; Blain 
et al. 1999; Blain et al. 2000; Eales et al. 1999, 2000; and 
Bargeretal. 1998, 1999). However, insufficient angular resolu- 
tion does not allow us to identify most of these objects in other 
wavelengths. Forthcoming projects such as SIRTF, ASTRO- 
F, SOFIA, Herschel Space Observatory, NGST, and ALMA 
would bring about revolutionary progress from this situation. 
Comparison of these data with theoretical models of infrared 
galaxies will be the key to understand the "dark side" of galaxy 
formation in the next decade of extragalactic astronomy. 

Theoretical modeling of galaxy formation and evolution has 
a long history. In optical/NIR wavelengths, most theoretical 
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modelings can be roughly divided into two categories: the 
so-called "backwards" approach and "ab initio" approach. In 
the former approach (Tinsley 1980; Yoshii & Takahara 1988; 
Fukugita et al. 1990; Rocca-Volmerange & Guiderdoni 1990; 
Yoshii & Peterson 1991, 1995; Pozzetti et al. 1996, 1998; 
Jimenez & KashUnsky 1999; Totani & Yoshii 2000), the lumi- 
nosity function of local galaxies is used as an input to normalize 
the number density of galaxies. The local properties of galax- 
ies such as multi-band colors and chemical properties are also 
used to construct a reasonable model of star formation history 
and luminosity evolution of galaxies based on the stellar pop- 
ulation synthesis method. The evolution is then probed back- 
wards into the past to predict observables such as galaxy counts 
and redshift distributions. The formation epoch and merging 
history of galaxies cannot be predicted in this framework, and 
hence they are introduced as phenomenological parameters that 
can be inferred from comparison with observational data. In the 
latter approach (Kauffmann et al. 1993; Cole et al. 1994, 2000; 
Somerville & Primack 1999; Nagashima et al. 2001), on the 
other hand, the formation epoch and merging history of galax- 
ies are predicted by the standard theory of structure formation 
in the CDM universe. In these models the local properties such 
as luminosity function are outputs of the model which should be 
compared with observations. However, although the formation 
and evolution of dark matter halos are rather well understood 
and can be properly predicted, our knowledge about baryonic 
processes such as star formation, supernova feedback, or galaxy 
merging is still very poor, and a number of phenomenological 
parameters must be incorporated, making the comparison of a 
model and observed data rather complicated. 

Both of these two approaches have been used also to make 
predictions of infrared galaxy counts and infrared CBR (see 
Franceschini et al. 1994 for the former, Guiderdoni et al. 1998 
and Devriendt & Guiderdoni 2000 for the latter, and Tan, Silk, 
& Balland 1999 for a somewhat hybrid approach between the 
two). However, because of our poor knowledge and under- 
standing of dust formation and emission, there are consider- 
able difficulties in theoretically predicting infrared luminosity, 
spectral energy distribution (SED), and their evolution of dust 
emission, compared with the optical/NIR modeling based on 
the stellar population synthesis method. This is the motivation 
of another "empirical" approach^ in the MIR and FIR wave- 
lengths. In this kind of approach, the infrared luminosity evolu- 
tion is introduced by some phenomenological functional forms 
and constraints on luminosity evolution is derived from com- 
parison with observed data (Beichman & Helou 1991; Oliver, 
Rowan-Robinson & Saunders 1992; Blain & Longair 1993, 
1996; Pearson & Rowan-Robinson 1996; Malkan & Stecker 
1998, 2001; Takeuchi et al. 1999, 2001a, 2001b; Roche & 
Bales 1999; Gispert, Lagache, & Puget 2000; Wang & Bier- 
mann 2000; Xu et al. 2001; Franceschini et al. 2001). 

There is yet another approach so-called "cosmic chemi- 
cal evolution", in which the universe is treated as a uniform 
medium and the mean gas consumption by stars in the whole 
universe is considered (Pei & Fall 1995; Fall, Chariot, & Pei 
1996; Pei, Fall, & Hauser 1999; Sadat, Guiderdoni, & Silk 
2001). The cosmic star formation history is determined by the 
redshift evolution of the cosmic mass density of neutral gas and 
the cosmic mean metallicity, which is inferred from quasar ab- 
sorption line systems. Although this is a simple and beautiful 



approach to predict CBR (but not galaxy counts), several prob- 
lems exist in this kind of approach, which could be serious espe- 
cially in modeling infrared radiation from dust. First, quasar ab- 
sorption systems are inevitably biased to gas-rich systems, i.e., 
systems in which star formation has been relatively inefficient, 
and hence their mean metalUcity does not necessarily represent 
the mean of the whole universe. For example, if an elliptical 
galaxy has formed at z ~ 4, this galaxy has been already in pas- 
sive evolution phase at z ~ 2. Then it can never be observed 
as a quasar absorption system because of complete exhaustion 
of interstellar gas, although significant stars and metals have 
been produced. Galaxies in dusty starburst phase may not be 
traced by absorption line systems, either, because background 
quasars may be completely extincted and cannot be observed. 
Therefore, the metallicity and star formation activity observed 
in absorption line systems must be an underestimate of the true 
mean in the universe. Furthermore, extinction by interstellar 
dust and reradiation in infrared bands are very sensitively de- 
pendent on dust opacity and its geometrical distribution within 
a galaxy, but this important information is completely missed. 

In this paper we construct a model of infrared galaxy counts 
and CBR based on the backwards approach for luminosity 
evolution of galaxies, by extending the model for optical/NIR 
galaxy counts (Totani & Yoshii 2000, hereafter TYOO) to in- 
clude dust emission component in infrared bands. After the 
work by Franceschini et al. (1994, F94 hereafter), only a lim- 
ited number of models have been pubUshed based on this ap- 
proach covering FIR bands, while there has been a dramatic in- 
crease of quality and quantity of observational data in the local 
universe as well as for faint/high-z galaxies. Although the back- 
wards approach has a disadvantage that formation epoch and 
merger history must be treated phenomenologically, it has an 
advantage that the number of model parameters is much fewer, 
and coimection of the model galaxies to the local galaxy pop- 
ulations is clearer, than the ab initio approach. Another reason 
for the relatively low activity of work in this direction is the 
difficulty of constructing realistic evolutionary models that can 
consistently describe both the optical and infrared data at the lo- 
cal universe as well as high-z, as argued by Franceschini et al. 
(2001). This is indeed what motivated more empirical studies 
confined to the infrared bands. 

The aim of this paper is to explain the bulk of optical and 
infrared data consistently by realistic evolutionary models of 
known and relatively normal galaxy populations at the local 
universe, without introducing hypothetical populations or para- 
metric description of luminosity evolution. An important im- 
provement of our model from the F94 model is that mass scale 
dependence is introduced for dust extinction and reradiation, 
based on the size-luminosity relation of galaxies observed in 
the optical bands. In the F94 model, there was no physical dif- 
ference in model galaxies of a given type with different masses; 
massive galaxies were simply a scaled-up version of smaller 
objects. We will show that this point is essential to understand 
various local properties of infrared galaxies. 

Another important characteristic of this work is modeling of 
the infrared SED evolution. In most of the previous infrared 
models, including ab initio and empirical approaches, the dust 
SED is empirically modeled in such a way that the dust SED is a 
one-parameter family specified by the total infrared luminosity 
of galaxies. The relation between the dust SED and luminosity. 
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which is observed in local infrared galaxies, is characterized by 
a gradual increase of characteristic dust temperature T^^st with 
infrared luminosity Ljr and is often expressed simply as a lin- 
ear relation between Tdust and logLiR (e.g. Smith et al. 1987; 
Soifer & Neugebauer 1991). However in this paper, we argue 
that this relation is not physically warranted for high-z galaxies. 
In a physical sense, the luminosity is extensive, i.e., a quantity 
which is dependent upon the amount of substance present in the 
system. On the other hand, the dust temperature is intensive, 
which is a specific characteristic of the system and is indepen- 
dent of the amount of material concerned. Therefore these two 
must not be related by only one sigle relation. At least there 
should be another extensive quantity such as the total mass of 
dust in a galaxy Mdust- In the empirical LiR-Fdust relation, the ex- 
tensive scale of an object is completely missed. The empirical 
J^iR-Tdust relation of local infrared galaxies should be considered 
as a projection of the plane between Ljr, Tdust, and Mdust- Here 
we will construct a model on physical basis to describe the in- 
frared luminosity as a two-parameter family of Mdust and Tdust- 
The dust temperature is physically calculated, without using the 
empirical LiR-Tdust relation. In addition to the thermal emission 
from big grain dust, the warm/nonequilibrium components of 
dust emission from small grain dusts and PAH features are also 
added to the model SED to be consistent with the observed in- 
frared SED of galaxies. 

It is important to check that the model can successfully re- 
produce the local properties of infrared galaxies, before pre- 
dicting the high-z quantities such as counts and CBR, since in 
our model all the inputs are the optical/NIR properties of lo- 
cal galaxies, and those in infrared bands are outputs predicted 
by the model. We will compare the model prediction with the 
observed properties of local infrared galaxies, such as the corre- 
lation between optical and infrared luminosity, the iiR-7dust re- 
lation, and infrared luminosity function. Indeed, we will show 
that our model gives a reasonable explanation for these rela- 
tions as well as the scatter along the mean relation. Recently 
Granato et al. (2000) have presented a very sophisticated semi- 
analytic model of infrared galaxy formation to calculate various 
properties of local infrared galaxies, although predictions for 
high-z galaxies have not been made yet. They have also found 
that most of the observed properties can be reproduced in their 
model framework. However, a large number of parameters have 
been introduced to treat various physical processes, and the ori- 
gin of the observed properties of local infrared galaxies has not 
been clearly discussed. Here we try to shed hght on the under- 
lying physics of these properties by a model with a much fewer 
number of adjustable parameters. 

Then we will make predictions of infrared galaxy counts and 
CBR by our backwards evolution model. We found that the em- 
pirical LiR-7dust relation breaks down at high-z, especially in the 
intense starbursts expected at the formation of elliptical galax- 
ies. The dust temperature should be much higher (40-80K) 
than predicted by mere extrapolation of the empirical LiR-7dust 
relation. We will discuss the implications of our results for the 
star formation history of the universe, especially the formation 
epoch of elliptical galaxies. Our result also has an interesting 
implication for the prima facie conflict (Stecker 2000) between 
the upper limit on MIR CBR from TeV gamma-ray observa- 
tions and the FIR CBR detections by the COBE satellite. The 
calculation of intergalactic optical depth of TeV garmna-rays 
based on our model will also be presented. 

The paper wiU be organized as follows. The model of galaxy 



evolution in the optical bands, on which the infrared modeling 
is based, will be summarized in §2. Then the extension of the 
optical/NIR model to include the dust emission is described in 
detail in §3. Our model predicts local properties of infrared 
galaxies, and they are compared with the observations in §4. 
The predictions of faint infrared counts and CBR will be pre- 
sented to be compared with observed data, in §5. Implications 
and discussions on our results are given for several topics in §6, 
and finally we summarize and conclude this paper in §7. Un- 
less otherwise stated, we will adopt a cosmological model with 
h = i/o/(100km/s/Mpc) = 0.7, = 0.2, I^a = 0.8. 

2. MODEL OF OPTICAL/NIR GALAXIES 

We use the same model as that used in TYOO, which is in 
overall agreement with the optical counts and photometric red- 
shift distributions in HDF (TYOO), and ^-band counts, colors, 
and size distributions in the SDF (Totani et al. 2001c). The 
detail of the model is given in TYOO and here we summarize 
important properties of the model. Although active galactic nu- 
clei (AGNs) may have some contributions to infrared galaxy 
counts and CBR, most of previous studies indicate that they are 
not the dominant population in counts or CBR (e.g. Rigopoulou 
et al. 2000; Barger et al. 2001). Therefore we will consider 
only normal galaxies without AGNs. Galaxies are classified 
into five morphological types of E/SO, Sab, Sbc, Scd, and Sdm, 
and their evolution of luminosity and spectral energy distribu- 
tions (SEDs) is described by a standard galaxy evolution model 
in which their star formation history is determined to repro- 
duce the present-day colors and chemical properties of galaxies 
(Arimoto & Yoshii 1987; Arimoto, Yoshii, & Takahara 1992). 
These models include chemical evolution giving gas fraction 
/gas and metallicity Z as a function of time and morphological 
type. We assume that the original dust-free SED of galaxies and 
its evolution do not depend on the galaxy luminosity (or mass). 
Their SED evolutions are shown in Fig. 1 of TYOO. The number 
density and morphological type mix are determined by the type- 
dependent present-day B-band luminosity function, and here 
we adopt that of SSRS2 survey (Marzke et al. 1998). Therefore, 
galaxies in our model are identified by morphological type and 
present-day B luminosity. All galaxies are assumed, for sim- 
plicity, to be formed at a single redshift, zf- We will adopt the 
formation redshift within a range of zf = 2-5 to see the sensi- 
tivity of our predictions on this parameter. 

A unique characteristic of the TYOO model is that the size 
and luminosity profile of galaxies are modeled and used to cal- 
culate the detection completeness of high-z galaxies in a realis- 
tic way based on the formulation developed by Yoshii (1993). 
Selection effects are very important in optical bands because 
galaxies are detected as extended sources and galaxy images 
are significantly affected by observational seeing and the cos- 
mological dimming of surface brightness. In infrared bands, 
angular resolution is generally much larger than typical appar- 
ent galaxy size, and hence they are essentially observed as point 
sources. Therefore, extended profile of galaxies is not relevant 
to the detection probability. However, the size of galaxies is 
still very important for the opacity of dust extinction. Here we 
summarize the modeling for galaxy sizes. We use a relation be- 
tween 5-band luminosity and effective radius to calculate a 
characteristic size of a galaxy as a function of its present-day 
B-band luminosity. The relation between the present B lumi- 
nosity Lb and effective radius r,, is assumed to be a power-law 
as re oc Lg^^'' for each galaxy type, and the normaUzation and 
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p are determined by fitting to the empirical relation observed 
for local galaxies. There is a considerable scatter in the empiri- 
cal r^-LB relation (about 0.22 in logr^), and this point wiU also 
be considered later. The observed LB-rg relations for elliptical 
and spiral galaxies and fitting parameters to them can be seen 
in Fig. 3 of TYOO. In this paper we assume that galaxy size 
evolves only when number evolution of galaxies takes place 
(see below). 

A simple picture for galaxy evolution is a so-called pure lu- 
minosity evolution (PLE) model in which there is no number 
evolution of galaxies to high redshifts. Although some number 
evolution is naturally expected in the hierarchically clustering 
universe dominated by cold dark matter, it is very difficult to 
construct a realistic model of merger history of galaxies espe- 
cially when one tries to predict infrared emission from dust. 
Therefore, we use the PLE model as our baseUne model. The 
difficulty of modeling realistic merger history of infrared galax- 
ies is uncertainty of size evolution of galaxies by the merger 
processes. In optical/NIR bands, size of galaxies is not so im- 
portant when extinction is not serious and the detection incom- 
pleteness is ignored. In infrared bands, however, emission from 
dust is always sensitive to dust optical depth that is directly cor- 
related with size of galaxies. Therefore size evolution induced 
by merging processes should have significant effect on the lu- 
minosity evolution of infrared galaxies. However, size evolu- 
tion in merging galaxies is hardly understood and very diffi- 
cult to construct a reliable model, even in the ab initio models 
where merging history of dark matter halos can be predicted 
rather well. Therefore, in this paper we will give only a simple 
calculation in the case of number evolution of galaxies, to grasp 
a rough behavior of the effects of number evolution. We intro- 
duce a simple merger model characterized by the evolution of 
the Schechter parameters of the luminosity function, in which 
total optical luminosity density is conserved, i.e., </>* oc (1 + 
and L%(x{\ +zT'^. For simplicity, we assume that all types of 
galaxies have the same number evolution. The size evolution is 
determined by assuming that stellar luminosity density within a 
galaxy is not changed in the merger process, i.e., a . How- 
ever, this is not warranted and should be considered as a simple 
toy model of merging galaxies. 

Absorptions by interstellar dust is calculated as follows. The 
dust opacity is assumed to be proportional to the gas column 
density and metalUcity, and it can be expressed by a function of 
restframe wavelength A, redshift z, morphological type denoted 
by k, and present-day B luminosity Lb^ as 

, 3,,,,, ( N^.s{z,k,LBfi) \ (Z{z,k)X 
Ti{X,z,k,LBfi)=X{X)Av,M^[ — , 

\ ^gas,MW / V -^MW / 

(1) 

where X{X) =A\/Ay is the extinction curve, A'gas the gas col- 
umn density, and Z the metallicity of interstellar gas. In this 
paper we use the Galactic extinction curve (e.g., Pei 1992), and 
s = \ as a standard prescription. On the other hand, depen- 
dence of the extinction curve on metallicity can be taken into 
account by power-law interpolations based on the solar neigh- 
borhood and Magellanic Clouds, with s = 1.35 for A < 2000A 
and s= 1.6 for A > 2000A (Guiderdoni & Rocca-Volmerange 
1987; Franceschini et al. 1994). We will also try this prescrip- 
tion later. Here the subscript MW denotes the value of the 
Milky Way at z = 0, and we assume that the Milky Way is an 
Sbc galaxy with Ly = 1.4 x lO'^L© (Binney & Tremaine 1987). 
The normalization of the optical depth is fixed by the parameter 
'^I'.MW, which will be determined later. The column density can 



be expressed as: 

yVgas(z, k, LBfi) (X Mb(k, LBfi) fg^siz, k) re(z, k, LsfiT^ , (2) 

where Mb is the total baryon mass of a galaxy and /gas the 
gas mass fraction. We use the present-day Lb^ and mass-to- 
luminosity ratio iMb/LBfi)k given by dust- free galaxy evolution 
models, to calculate Mb. In a strict sense, using {Mb/Lsfi) 
of a dust-free evolution model is an approximation, because 
the observed o is already obscured by dust. However, we 
use this approximation for the simplicity of numerical calcula- 
tions. This approximation is reasonable provided that the opti- 
cal depth of extinction at B-band in present-day galaxies is not 
much greater than unity, which is supported by observations 
for the majority of nearby galaxies. However, as we will argue, 
more massive or luminous galaxies are more dusty and hence 
the above estimate of gas column density may be underestimate 
for the brightest galaxies. This point will be discussed later. 

Extinction of galactic light heavily depends on the spa- 
tial distribution of dust in a galaxy. Here we adopt two ex- 
treme models of spatial distribution of dust: the intervening 
screen model and the slab (i.e., the same distribution for stars 
and dust) model. The attenuation factor of galactic light is 
f-MiX.z.k.LBfl) = exp(-rrf) and {1 - exp(-Td)} / for the for- 
mer and the latter, respectively. Then we obtain the extincted 
SED of a model galaxy as 

LxiX,z,k,LBfi) = f^f^iXB,0,k, LBfi) fMiX,z,k, LBfi) 

X Ll'%X,z,k,LBfi) , (3) 

where L^^^ is the SED of dust-free model galaxies, and the first 
factor of accounts for the fact that the observed Lba is al- 
ready extincted. The difference between the screen and slab 
models is not significant for galaxies that are not so dusty (i.e., 
optical depth < 1), but is quite significant for dusty galaxies. 
Although the slab model seems to be more reasonable, several 
observations of starburst galaxies cannot be explained by the 
slab model. The observed correlation between power-law index 
of UV spectra and Balmer line ratio of starburst galaxies indi- 
cates that the observed reddening of starburst galaxies is larger 
than expected from the slab model, and at least some fraction of 
dust should behave like a screen (Calzetti, Kinney, & Storchi- 
Bergmann 1994, see also Gordon, Calzetti, & Witt 1997). Al- 
though the screen model may seem unreasonable, it is rather 
reasonable that a fraction of dust behaves like a screen with a 
distribution more extended than stars, since strong wind from 
starbursts could blow out the dust particles along with interstel- 
lar gas. This could occur in the scale of star forming regions 
within a galaxy and/or the whole scale of a starbursting galaxy. 
Some ultra luminous infrared galaxies (ULIRGs) are known to 
have infrared luminosity which is more than 100 times higher 
than optical/NIR luminosity. In the slab model, this requires an 
enormous amount of dust with > 100, while in the screen 
model dust amount can be more modest. Therefore the screen 
model may be better than the slab model as a phenomenolog- 
ical prescription. We then use the screen model as a standard, 
and also use the slab model to see the dependence on the dust 
distributions. 

The model described here has already been compared com- 
prehensively with the counts and redshift distributions of galax- 
ies observed in optical and NIR bands (Totani & Yoshii 2000; 
Totani et al. 2001c). It has been found that all data can be rea- 
sonably explained by this model assuming a A-dominated flat 
cosmology. The NIR data can be explained weU by the PLE 
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model (Totani et al. 2001c), while the optical data suggest some 
modest number evolution (r? ~ 1) (Totani & Yoshii 2000). The 
apparently discrepant results between optical and NIR bands 
may reflect dependence of number evolution on galaxy types, 
i.e., almost no number evolution for elliptical galaxies while a 
stronger number evolution for late-type galaxies (see Totani et 
al. 2001c for more detailed discussion). 

3. EXTENSION TO INFRARED GALAXIES 

3.1. Determination of Effective Dust Temperature 

In the optical modeling described above, we have already 
calculated the absorption of stellar light by dust in a galaxy of 
given type at any redshift (or time). Therefore we can straight- 
forwardly obtain the total infrared luminosity emitted by dust 
particles as 

j dXf;^iXB,0,k,LB,o) 

X [l-fA\z,k,LB,o)]Ll''(X,z,k,LB,o). (4) 

Then the task remained to predict infrared flux at a given ob- 
serving band is to model the SED of dust emission. In a number 
of previous papers, an empirical relation between the infrared 
SED (or characteristic dust temperature Tdust) and Lir has been 
utilized here. However, as mentioned in Introduction, it is phys- 
ically unreasonable to assume that the luminosity is determined 
by only one quantity, i.e., temperature. The total luminosity 
is an extensive quantity that is proportional to dust mass, pro- 
vided that galaxies are transparent for dust emission. If a galaxy 
is simply scaled up in mass, the infrared luminosity should also 
increase but the temperature should not, which is an intensive 
quantity determined by the local physical state of dust particles. 
We have no reasonable physical explanation for this empirical 
relation so far, and it is not warranted at aU that this can be 
appUed to high-redshift galaxies. 

Rather, dust temperature should be determined by an energy 
balance condition; i.e., temperature should be a value at which 
dust emissivity is consistent with total amount of stellar light 
absorbed by dust. We assume that total infrared luminosity is 
dominated by thermal radiation of big grain dust, which is well 
described by the so-called graybody spectrum, i.e., blackbody 
spectrum multiplied by the dust emissivity that is described by 
a power of frequency as e oc i^''', in the FIR regime. This is a 
good approximation for local galaxies, and we will find that it 
holds also for high-z galaxies whose dust temperature is higher. 
Then the energy balance condition can be written as: 



Tdu 



4+7 



(5) 



where Mjust is the total dust mass in a galaxy and the dust emis- 
sivity index 7 is typically 7 = 1-2 (Boulanger et al. 1996; Hirao 
et al. 1996; Calzetti et al. 2000; Dunne et al. 2000). (The case 
of 7 = corresponds to the blackbody spectrum.) Throughout 
this paper we use 7 = 1.5. The factor /d is a mass fraction of 
dust heated by stellar light. In the case of the slab dust, all dust 
should be heated equally, while in the case of the screen dust 
only a part of dust can be heated by penetrating fight. There- 
fore, 



/d 



1, (slab) 

l-exp[-rd(Aeff)] 



Td(Aeff) 



-, (screen) 



(6) 



where the characteristic wavelength Aeff of absorbed stellar 
light is defined as: 

JdXX[l-fMX)]Lf%X) 



A. 



(7) 



/^/A[l-/a„(A)]Lfr-(A) 
The subscript * denotes the reference point to fix the normal- 
ization of eq. (5). We take this reference point as Lir^ = 
lO^°h~^LQ and Jdust,* = 19 K for a galaxy type Sbc at z = 0, to 
be consistent with the observed Ljr and Tdust relation. The rel- 
ative dust mass ratio (Mdust/A^dust,*) can be calculated for any 
model galaxy assuming that dust-to-metal ratio is constant, i.e., 
^dust oc /gasMbZ. In this way we obtain Tdust for any model 
galaxy. 

3.2. The Infrared SED Templates 

The overall shape of a galaxy SED is basically determined 
by dust size spectrum and strength of radiation field. Here we 
should note that the entire range of infrared SED of dust emis- 
sion cannot be described by a simple graybody spectrum with a 
single temperature. While FIR peak of the infrared SED is well 
described by a single-temperature graybody spectrum of cold, 
big grain dust, the MIR emission is dominated by radiation 
from heated small grains and PAHs that are not in thermal equi- 
librium with ambient radiation field. Because of their small heat 
capacity compared with the heating photon energy, these dust 
grains are heated stochastically by one or two photons, and their 
'temperature' temporally variates violently (e.g., Purcell 1976; 
Draine & Anderson 1985; Siebenmorgen, Kriigel, & Mathis 
1992; Draine & Li 2001; Li & Draine 2001). This phenomenon 
makes the MIR continuum very broad, and causes a power-law 
Uke spectrum. Thus, it may be improper to describe the MIR 
continuum as a superposition of multi-temperature blackbody 
spectra. 

Though the MIR continuum is dominated by the contribu- 
tion from non-equilibrium dust grains, it is well correlated with 
temperature of big grain dust, Tdust (Shibai, Okumura, & Onaka 
2000; Dale et al. 2001). It enables us to construct the over- 
all MIR-FIR SEDs as a unique function of Tdust- We utilize 
the observed correlations between MIR band fluxes (Dale et al. 
2001) and construct the MIR regime of the model SED in an 
empirical way. This is a standard prescription used for the dust 
SED modeling in a number of previous papers. This is equiva- 
lent to assuming that the local physical state of dust particles is 
specified by one parameter, i.e., characteristic temperature. It 
is physically not unreasonable, though not warranted, to expect 
that this holds also for high-z galaxies, unless the composition 
or size distribution of dust particles in high-z galaxies is dras- 
ticaUy different from local galaxies. The point different from 
previous models is that we do not simply relate Tdust to the in- 
frared luminosity Lir, but derive Tdust by the energy balance 
condition. 

Details of the construction of the MIR SED are as follows. 
The spectrum of thermal, big grain dust particles is described 
by a graybody spectrum with Tdust. The MIR continuum com- 
ing from small particles is described by a broken power law 
with two different indices, which are determined by the ob- 
served fiux at 12 and 25 /um compared with that in 60 fim at 
which the graybody is dominant. These relations are: 

log ^ = -0.026 Tdust+0.367 , (8) 

-J 60 

and 



log^=-1.751og 1^-0.866. 

^25 'J 100 



(9) 
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where is the flux density of infrared galaxies in units of 
Wm"^Hz~' at A = X ij.m, and the temperature Tdust is in unit 
of K. In addition to this, the PAH features are taken from spec- 
troscopic data obtained by ISO SWS (Dale et al. 2001), and 
their strength is delcrmined by the flux at 6.75 /um compared 
with the continuum emission as: 



log 



S6.- 



-1.571og56o/5ioo-0.618 (Seo/Sm > 0.4) 
-0.104 



(10) 



5i5 -0.104 (W5ioo<0.4) 

Thus we obtain the MIR-FIR SEDs of galaxies with various 
dust temperature. Here we sunnmarize the shape and frequency 
range of each component as follows: i) a power-law compo- 
nent (1.0 x lO'^-l.Ox 10'^ [Hz]). The longer wavelength edge 
of this component is set to decline as a power law oc v^^'', and 
shorter wavelength edge is set to decline exponentially, ii) a 
PAH continuum obeying a power law (1.0 x 10^^-3.0 x lO" 
[Hz]). This component is connected to the above one smoothly, 
iii) a PAH band emission component (2.5 x lO'-'-l.O x 10'"^ 
[Hz]). Finally we added radio continuum to the IR SED 
along with the tight radio-FIR correlation (e.g. Helou, Soifer, 
& Rowan-Robinson 1985; Condon 1992), 

log5i.4GHz = log(0.86756o + 0.3365ioo)-2.3 . (11) 

Figure 1 shows our model SED for various values of T^ast- A 
trend is seen that the MIR component coming from nonequilib- 
rium dust emission becomes weaker compared with the gray- 
body emission from big grain dust, with increasing tempera- 
ture. This trend is reasonable because the high temperature is 
equivalent to stronger flux of heating radiation field. When the 
heating intensity is stronger, the heating photons are more fre- 
quently absorbed by small dust grains, and hence more dust 
particles become in thermal equihbrium, making the spectrum 
closer to the graybody. There is also a large body of evidence 
that the PAH are likely to be destroyed, or that their emission 
is severely damped, in regions of high heating intensity (Dale 
et al. 2001 and references therein). These trends have an im- 
portant implication for the compatibiUty of MIR CBR and TeV 
gamma-ray observations, which will be discussed in §6.4. 

Figure 2 shows the model infrared SEDs for T^ust = 21 and 
45 K, corresponding to the cases of the Milky Way and M82. 
Our model is in good agreement with the observed SEDs in 
both cases. Figure 3 is the same, but showing the fit of our 
model to other four infrared galaxies with various values of 
Tdust- These results indicate that our model of infrared SEDs 
is in reasonable agreement with observed infrared galaxies in a 
wide range of dust temperature. 

Now we can predict various observable quantities of infrared 
galaxies, i.e., local properties as well as faint galaxy counts and 
CBR composed by high-redshift galaxies, based on the model 
of optical galaxies described in the previous section. The only 
one parameter not yet determined is Ay_MW, the overall normal- 
ization of dust extinction. Increasing Av.mw makes galaxies 
more dusty and hence more luminous at infrared wavebands. 
Infrared galaxy counts at a fixed flux then increase with Avmw- 
In this paper we set this parameter to a value at which bright 
infrared galaxy counts at 15, 60, and 90 /zm are consistent with 
the IRAS counts (see §5.2). This is found to be Av^mw ^ 0.7 
for our baseline model. However, it is important to check 
whether the model is consistent with observations for local in- 
frared galaxies, before we discuss the result of counts and the 
CBR. This wiU be presented in the next section. 

4. LOCAL PROPERTIES IN INFRARED: MODEL VERSUS 
OBSERVATIONS 



Here we compare our model predictions with the proper- 
ties of local galaxies observed by IRAS. We use the sample 
of the PSCz survey of 1541 1 IRAS galaxies selected by 60 fjm 
flux (Saunders et al. 2000). Out of this sample, we further se- 
lected 12324 galaxies detected in the 100 //,m band as well, for 
which the bolometric infrared luminosity can be rehably esti- 
mated. Systematic and reliable optical flux measurement is not 
available for all of these galaxies (see the ReadMe file of the 
PSCz catalog of Saunders et al. 2000 for detail). We use flie 
Zwicky/RC3 magnitude as B-band flux, which is available for 
5803 optically bright galaxies. 

4.1. Optical/Infrared Luminosity Ratio 

Figure 4 shows an Lg iyFi^ at 4400A) versus Lir (bolomet- 
ric) plot. The data points are the 5803 galaxies with available 
B-band magnitudes. A clear trend can be seen that more lumi- 
nous galaxies are more dusty, i.e., showing high Lir /Lb ratio. 
On the other hand, the thick solid and dashed lines are our pre- 
diction averaged over all the four spiral galaxy types at z = 0, 
for the screen and slab dust, respectively. Here galaxy types are 
averaged by their relative proportions calculated from the type- 
dependent B-band luminosity function. It should be noted that 
elliptical galaxies are not included in this prediction because 
they do not have dust and hence infrared emission at z = in our 
model. The trend of increasing (Ljr/Lb) ratio with luminosity 
is well reproduced by the model. This trend stems from the 
adopted B luminosity versus size relation; dust opacity, which 
is proportional to the metal column density, increases with Lb 
when the locally observed LB-rg relation is taken into account. 
Consequently, more massive galaxies emit their radiation more 
in infrared bands. There is a considerable scatter in the Lg-r^ 
relation, and this should produce scatter in the optical depth and 
hence a scatter along the mean Lir-L^ relation. The thin-sohd 
and thin-dashed curves are the predictions when the L^-re rela- 
tion is shifted by ±lcr = 0.223 in log r^. from the mean relation, 
where the line markings are the same with those of thick lines. 
The scatter is also in reasonable agreement with the observed 
scatter. 

There is some offset between the prediction and mean ob- 
served LiR-Lfi relation. We suggest that it is due to the incom- 
pleteness of the sample; the RC3 B-magnitude is available only 
for galaxies bright in optical bands, which are about half of the 
whole infrared sample of 12324 galaxies. The galaxies without 
B-magnitudes are then expected to be less luminous in optical 
bands. This effect could lead to the lack of infrared luminous 
galaxies compared with the model prediction in Fig. 4. 

4.2. Dust Temperature versus Infrared Luminosity 

Figure 5 shows the infrared bolometric luminosity versus 
dust temperature plot. The data points are the 12324 galaxies in 
the IRAS PSCz sample with 100-/xm flux available. Dust tem- 
peratures of the JRAS galaxies are estimated from their Seo/Sm 
flux ratio based on the temperature measurement of the Galactic 
dust by IRAS and COBE (Nagata et al. 2001, in preparation). 
The thick/thin solid and dashed lines are the model predictions, 
with the same line markings corresponding to model predic- 
tions shown in Fig. 4. A trend of gradual increase of mean T^ast 
with Lir, as well as the scatter along the mean relation, is again 
well reproduced by the present model. This also stems from 
the fact that more luminous galaxies are more dusty, as argued 
in the previous subsection. This fact requires that Lir should 
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increase faster than a simple scaling to Mjust, when the mass 
of galaxies is increased. Considering the energy balance of 
LiR oc Mdust7dus7' the dust temperature must gradually increase 
with Mdust or L,r. 

There is a clear lack of low temperature galaxies with Tjust ^ 
16K in the IRAS sample. We suggest that this sudden lack 
may come from selection effects. As we mentioned above, the 
dust temperature is calculated based on the 60-/im and 100- 
fim fluxes. However, the sample is selected at IRAS 60-/xm 
bandpass, and very cold galaxies have an IR SED peak at wave- 
length much longer than 60-/im. In addition, such galaxies have 
low luminosity in general. Therefore, cold galaxies are easily 
dropped by the sampling procedure in the sample used here. 
This is likely to be the cause of the sudden cutoff of the data on 
the diagram. 

4.3. Infrared Luminosity Function 

Finally we compare the model prediction of local infrared 
luminosity function (at 60/xm) to the estimated one from ob- 
servations (Saunders et al. 1990), in Fig. 6. It is well known 
that the infrared luminosity function is considerably different 
from that of optical galaxies that is well fitted by the Schechter 
function. The dot-dashed line is the Schechter-type luminosity 
function translated from the fi-band luminosity function assum- 
ing vLv{B)/uL,^{6Qijm) = 1. In sharp contrast to the rapid ex- 
ponential decline of the optical luminosity function above L*, 
infrared luminosity function extends well beyond Lgo ~ lO'^L© 
where Leo is vL^ at 60/(xm (Fig. 6). Despite this intrinsic dif- 
ference between the two, our model prediction with the screen 
dust, which is based on the B-band luminosity function having 
the Schechter form, well reproduces the shape of 60/im lumi- 
nosity function up to Lgo ~ 10' ^Lq. Here, the infrared luminos- 
ity function is calculated as 

i-f, V M^B), (12) 

where (j)B,k is the S-band luminosity function of the k-ih type 
galaxies. This is again due to the fact that more luminous galax- 
ies are more dusty; when optical luminosity is increased, in- 
frared luminosity increases more rapidly than that, and hence 
an optical luminosity range is significantly broadened when it 
is translated into an infrared luminosity range. This effect is fur- 
ther strengthened for the infrared luminosity function at 60/im, 
since the gradual increase of temperature with bolometric lumi- 
nosity makes the increase of 60/im luminosity even faster. 

The models fail to reproduce the data at ultra-luminous re- 
gion of L > lO'^L©, possibly because of the contribution of pe- 
cuhar starbursting galaxies under galaxy interactions, or AGNs. 
In fact observations suggest that the fraction of AGNs increases 
with infrared luminosities (Veilleux, Kim, & Sanders 1999). 
Our model does not take into account these rare populations of 
ULlRGs at the local universe, since their contribution to counts 
or CBR is expected to be small (Rigopoulou et al. 2000; Barger 
et al. 2001). Instead, we try to explain bulk of the infrared data 
by the population of normal galaxies at the local universe and 
their ancestors at high redshifts. It should also be noted that the 
approximation about mass-to-hght ratio (Mb /Lb) breaks down 
for very bright and dusty galaxies, as mentioned in §2. There- 
fore our calculation is quantitatively unreliable and could be an 
underestimate at the brightest part of the infrared luminosity 
function. 

Recently the dust mass function of local galaxies is derived 
by Durme et al. (2000), based on infrared- submilUmeter SEDs. 



They found that the dust mass function can be fitted by the 
Schechter function, in contrast to the 60/im luminosity func- 
tion. This is also consistent with the expectation of our model; 
the dust mass is roughly proportional to the B-band luminos- 
ity unless the optical light is heavily extincted, and hence the 
dust mass function should be the same shape with the optical 
luminosity function. 

5. HIGH REDSHIFT INFRARED GALAXIES 

5.1. Evolution of Infrared Luminosity and SED 

In the following part of this paper we will make predictions 
for high-z infrared galaxies such as faint galaxy counts and 
CBR. First we show the evolution of infrared luminosity and 
big-grain dust temperature (Tjust) of model galaxies as a func- 
tion of time after the formation in Fig. 7. These models are con- 
structed to reproduce the observed optical SEDs and chemical 
properties at the typical age of present-day galaxies (>10Gyr). 
For Tdust the three cases of dust temperature modeling are shown 
in the three panels: the screen dust, slab dust, and using the em- 
pirical relation between Lir and Tdust of the local galaxies. Here, 
we used the best-fit relation to the observed data as: 

l.llogio(W'L0) + 8.5K, (13) 

for the empirical LiR-Tdust relation (shown as a dotted Une in 
Fig. 5). It should be noted that our evolution model depends 
on the mass scale of galaxies. Here we are showing the models 
for galaxies whose present-day 5-band luminosity is Mb = -20, 
as typical giant galaxies. It can be seen that dust temperature 
evolution is much stronger than expected from the empirical 
LiR-Tjust relation, especially for primordial elliptical galaxies. 
The peak dust temperature during the starburst phase of ellipti- 
cal galaxies in the slab-type dust model is ~30^0K, and even 
higher temperature of ~40-80K is found for the screen dust 
model, while the the empirical LiR-7dust relation predicts T^ust ^ 
25K in all time. This must have significant effect on the infrared 
SED of primordial elliptical galaxies, and hence on faint counts 
and CBR. 

In the primordial elUptical galaxies, the dust opacity is so 
high and hence Td ^ 1 for most of optical light. Therefore, 
the total absorbed light, i.e., Ljr is not much different between 
the slab and screen models. Rather, the difference of the two 
mainly comes from the difference of amount of dust irradiated 
by stellar radiation. In the case of sceen dust model with rj » 1 , 
only /d ^ Td(Aeff)~' of total amount of dust can be heated, while 
all dust is heated in the slab model (eq. 6). Consequently the 
temperature of the screen dust model must be higher to achieve 
higher luminosity per unit dust mass. It should be noted that 
the temperature estimate in the screen model is quite robust, in 
a sense that it does not depend on the total amount of dust. Once 
Td(Aeff) ^ 1 is achieved, the dust-mass dependence is canceled 
between /d and Mdust in the energy balance equation (eq. 5). 
Later, we will find that the screen model shows much better 
agreement with counts and CBR, than the slab model. 

The high temperature which we obtained here for primor- 
dial elliptical galaxies should be compared with some observa- 
tional estimates of dust temperature for high-z as well as local 
starburst galaxies. It should be noted that the dust tempera- 
ture obtained from a fit to observed spectra depends on several 
other parameters, such as the emissivity index 7 and opacity of 
a galaxy to FIR radiation. A same observed SED may be fit- 
ted with a higher Tdust when a smaller value of 7 is assumed 
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or the opacity to infrared radiation (tir) is significant. On the 
other hand, the model temperature derived here is an effective 
value derived from the total energy balance between dust mass 
and bolometric luminosity assuming 7=1.5 and low opacity to 
infrared light. Therefore, our result should be compared with 
observational estimates of temperature derived under the same 
assumptions of 7 ~ 1.5 and low opacity to FIR light. Keep- 
ing this point in mind, known high-z luminous starburst galax- 
ies seem to have a mean temperature of about 5 OK, under the 
above conditions for 7 and infrared opacity, and some of them 
show a high temperature of ~ 70K (Ivison et al. 1998a, b; Ben- 
ford et al. 1999). Recently a comprehensive study for the SED 
of local ULIRGs is published (Klaas et al. 2001), and it again 
shows that the ULIRG dust temperature ranges from ~40 to 
~80K, with 7 ^ 1.6 and TiR(lOOAim) - 0.5-5. These obser- 
vations then indicate that the temperature derived by our model 
for primordial elliptical galaxies is in a reasonable range of typ- 
ical known starburst galaxies, although it is much higher than 
the simple extrapolation of the mean LiR-Tjust relation of local 
galaxies. 

5.2. Infrared Galaxy Counts 

Figure 8 shows predictions of infrared galaxy counts com- 
pared with available data at 15, 60, 90, 170, 450, and 850 ^m 
wavebands by our basehne model. In the baseline model, we 
assumed a A-dominated flat cosmology with {h, Oq, ^^a) = (0.7, 
0.2, 0.8), formation redshift of z/? = 3 for all galaxy types, pure 
luminosity evolution without number evolution, and screen dust 
model. See figure caption for the references of observed data. 
The parameter of overall normaUzation of dust optical depth, 
A-v,MSN is determined here so that the count predictions fit best 
the brightest 7i?AS counts at 15, 60, and 90 /^m bands. This nor- 
malization has already been used in the previous sections. Our 
model predicts lower and higher counts than the bright counts 
observed in 15 and 60 /zm, respectively, and it may suggest 
that our modeling of infrared SED might be imperfect, but it 
might also reflect some systematic uncertainties in the calibra- 
tion of the observed flux in these wavebands (see e.g., Mazzei 
et al. 2001; Lari et al. 2001; Serjeant et al. 2001).' This fig- 
ure shows contributions from each galaxy type as well as the 
total counts. Nearby spiral galaxies are dominant in the bright- 
est counts, while contribution from elliptical galaxies rapidly 
appears when the flux is decreased down to a certain value, be- 
cause of the strong dusty starburst of these galaxies at their for- 
mation stage. Appearance of high-z starburst galaxies revealed 
by SCUBA in submillimeter wavelengths is in good agreement 
with count predictions of forming elliptical galaxies at z ~ 3 
in our model, and hence it can be interpreted as the emergence 
of forming elliptical galaxies. The redshift of z ~ 3 for these 
galaxies is also consistent with redshift estimates of SCUBA 
sources (e.g., Dunlop 2001). On the other hand, our model sug- 
gests that relatively weak evolution in the 15/xm-band counts 
at 5 < 10"^ Jy is not due to elUptical galaxies, but induced by 
evolution of spiral galaxies. Strong evolution suggested by ISO 
observations (e.g., Efstathiou et al. 2000; Matsuhara et al. 2000) 
in the 90/xm band may not be explained by the emergence of el- 
Uptical galaxies in our model. However, contribution from such 
galaxies with smaller zf, or bulges of spiral galaxies, could ex- 
plain the observed evolution (see Fig. 10 below). (It should be 
noted that the bulge and disk of spiral galaxies are modeled as 

' Calibration in the MIR and ¥]R requires state-of-the-art techniques and especial 
remains a significant uncertainty. We do not go further into the problem here. 



a mixture in our model.) 

Figures 9, 10, and 11 show the count prediction for the to- 
tal of aU galaxy types in cases of different cosmological mod- 
els, formation redshifts, and number evolutions, respectively. 
The dependence on cosmological models becomes important at 
flux much fainter than the present sensitivity limits, except for 
15/im. Therefore the uncertainty in cosmological parameters 
does not affect our conclusions significantly. Although there 
may be some uncertainties in our modeling of number evolu- 
tion, it seems that the counts of SCUBA sources in submillime- 
ter bands are better explained by primordial elhptical galaxies 
without number evolution, rather than the cases of number evo- 
lution. This may indicate that elliptical galaxies have formed 
by a single starburst at high redshifts, rather than by many star- 
bursts of smaller building blocks which would eventually merge 
into a giant elUptical galaxy after the starbursts. This is consis- 
tent with a number of literature about observational constraints 
on number evolution of elliptical galaxies (e.g., Totani & Yoshii 
1998; Im et al. 1999; Schade et al. 1999; Benitez et al. 1999; 
Broadhurst & Bouwens 2000; Daddi, Cimatti, & Renzini 2000; 
Totani etal. 2001c). 

Figure 12 shows some model predictions with different 
model prescriptions: using empirical LiR-7dust relation that has 
been used in most of previous studies (dotted Une), using slab- 
type dust geometry rather than the screen dust (dashed), and 
using a different prescription for dust opacity by Guiderdoni & 
Rocca-Volmerange (1987) as discussed in §2, rather than as- 
suming a constant dust-metal ratio (dot-dashed line). The last 
prescription gives almost the same prediction with our base- 
line model, and hence difference of the two prescriptions for 
dust/metal ratio is not important. On the other hand, the empir- 
ical LiR-7dust relation gives quite different prediction especially 
for elliptical galaxies at the dusty starburst phase, because the 
empirical relation gives considerably lower temperature than 
our baseline model as shown in the previous section. Conse- 
quently, counts at 90 and 170 /im are decreased, while those 
at 450 and 850 jim are enhanced. The slab dust has the same 
effect, although it is weaker than the case of the empirical Lir- 
Tdust relation. 

5.3. Cosmic Background Radiation 

We show our prediction of the CBR in optical and infrared 
wavelengths in Fig. 13, by our baseline model whose counts 
have been presented in Fig. 8. The reported detections of CBR, 
lower limits coming from integration of galaxy counts, and up- 
per limits on MIR CBR from TeV gamma-ray observations are 
also shown (see figure caption for references). The model is in 
overaU agreement with all available data. 

In optical/NIR bands, the faint-end slope of galaxy counts is 
very flat and hence the integration of galaxy counts with extrap- 
olation down to even fainter magnitude is convergent. There- 
fore the integration of galaxy counts in these bands should not 
be much different from the true CBR from galaxies. Since our 
model is based on a model of galaxy counts in optical/NIR 
bands which is in agreement with observations, it is trivial 
that our model prediction of optical/NIR CBR is in agreement 
with observational estimates based on integrations of observed 
galaxy counts. Such estimates are shown by error bars without 
symbols in Fig. 13, in which the contribution to CBR by galax- 
ies missed in deep surveys by various selection effects is taken 

scrutiny, and although a lot of efforts by experts have been devoted to it, there stiU 
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into account (Totani et al. 2001a). However, our model under- 
predicts the UV CBR flux compared with the observed count 
integrations, and this suggests that the dust extinction is too 
large in the model. In fact, in this paper we have set the overall 
normalization of dust, Avmw ^ 0.7, so that the bright infrared 
galaxy counts of our model are consistent with observed data. 
This is considerably larger than the typical Galactic extinction 
in relatively high latitude, Avmw ~ 0.2, which was used in our 
previous papers on optical/NIR counts (Totani et al. 2001a, b, 
c). This discrepancy is presumably coming from the scatter of 
dust optical depth in galaxies. In our model we apphed a sin- 
gle relation between B luminosity and galaxy size for the sim- 
plicity, but in reality the scatter along the mean L^-r^ relation 
should produce significant scatter of extinction, as discussed in 
§2 and §4. The infrared counts are expected to be dominated by 
galaxies with relatively large extinctions, while the optical/NIR 
counts, especially U band, would be dominated by those with 
less extinctions. This is probably the reason our prediction of 
UV CBR is lower than the observed count integration. Since 
the primary interest of this paper is in infrared bands, we will 
ignore the discrepancy in UV CBR in the foUowing part of the 
paper. 

There are several reports of detection of diffuse NIR CBR 
flux in the J (1.25/xm), K (2.2/xm) and L (3.5/Ltm) bands as 
shown by diamonds in Fig. 13 (Gorjian et al. 2000; Wright 
& Reese 2000; Wright et al. 2001; Matsumoto et al. 2001; 
Cambresy et al. 2001). All fliese reported diffuse NIR CBR 
fluxes are systematically higher by a factor of several than the 
integration of galaxy counts in the same bands. Our model is 
made to reproduce the observed galaxy counts, and hence our 
CBR flux agrees with the galaxy count integrations, and does 
not agree with the reported diffuse CBR flux. The contribu- 
tion from galaxies missed in deep surveys by selection effects is 
very unlikely to reconcile this discrepancy (Totani et al. 2001a). 
It may suggest some systematic uncertainties in diffuse CBR 
measurements, or the contribution to the NIR CBR from yet 
unknown sources, but it is beyond the scope of this paper. On 
the other hand, there is a recent report of detection of diffuse 
CBR in optical bands (Bernstein, Freedman, & Madore 2002). 
Although the reported diffuse flux again seems higher than the 
count integration estimations by Totani et al. (2001a), they are 
consistent with each other at la when all error bars (statistical 
and systematic of the diffuse measurements, and contribution 
from faint galaxies missed by deep surveys) are taken into ac- 
count. 

Figures 14, 15, and 16 show the CBR prediction for the to- 
tal of all galaxy types in cases of different cosmological mod- 
els, formation redshifts, and number evolutions, respectively. 
These figures show that the dependence on cosmological mod- 
els or number evolution is not large compared with observa- 
tional constraints, and it does not affect our conclusions sig- 
nificantly. On the other hand, the formation redshift can be 
constrained better since the sub-mm CBR flux rapidly increase 
with increasing zf due to the redshift effect on the FIR peak of 
dust emission SED, while there is an observational constraint 
by COBE/FIRAS. Our model suggests that the major formation 
epoch of elUptical galaxies should not be larger than zf ^ 5. 

Figure 17 shows some model predictions with different 
model prescriptions, with the same models and line markings 
with those in Fig. 12. There is again a remarkable differ- 
ence between the baseline model and the model using the em- 
pirical Z-iR-Tdust relation; the latter model seriously overpre- 



dicts the sub-mm CBR flux compared with the detection of the 
COBE/FIRAS. When the empirical LiR-7dust relation is used, 
temperature of starbursting elUptical galaxies does not become 
as high as our baseline model, and hence the redshift effect sig- 
nificantly pushes FIR peak of CBR toward longer wavelengths. 
In fact, because of this reason, most of the previous papers pub- 
lished recently have claimed that too intense starbursts weU be- 
yond z ~ 1 are disfavored by the FIRAS data. However, our 
model suggests, based on a physical argument, that high-z star- 
bursts of forming elliptical galaxies should have much higher 
temperature than previously expected from the empirical Lir- 
Tdust relation, and hence formation of elliptical galaxies at high 
redshifts (zp ^ 3) with dusty starbursts is allowed, or even fa- 
vored, from the data. This is consistent with the standard pic- 
ture of elliptical galaxy formation (Larson 1974; Arimoto & 
Yoshii 1987), i.e., intense initial starburst at high redshift fol- 
lowed by passive evolution to the present. The screen dust 
model seems much better than the slab-type model, as shown in 
Fig. 12 and 17. This has an interesting implication for the op- 
tical/NIR colors of primordial elliptical galaxies in dusty star- 
burst phase, and we wiU discuss it in the next section. 

6. IMPLICATIONS AND DISCUSSION 
6.1. Forming elliptical galaxies in optical/NIR bands: HEROs 

We have shown that the screen-like dust extinction gives 
quantitatively much better results both for the infrared counts 
and CBR than the slab-hke dust, if elliptical galaxies formed 
by starbursts at high redshifts of Zf > 2 as generally believed. 
This is consistent with the properties of nearby starburst galax- 
ies showing strong reddening which cannot be explained only 
by the slab-like dust. It also seems reasonable that at least a 
part of dust particles is blown out by strong galactic wind ex- 
pected in the formation stage of elliptical galaxies, to have a 
distribution similar to a screen. An implication of this result 
is that forming elUptical galaxies should show very red colors 
when observed in optical/NIR wavelengths. When dust opac- 
ity becomes much greater than unity, colors of galaxies could 
infinitely increase by screen dust, while colors should reach an 
asymptote in the case of slab dust since emission from stars 
located at the outermost region becomes dominant. 

In fact, such objects having unusually NIR colors oiJ—K> 
?)-A were discovered recently in the Subaru Deep Field (Mai- 
haraetal. 2001; Totani etal. 2001b). These colors are even red- 
der than those of typical extremely red objects (EROs) selected 
by optical-NIR colors, and cannot be explained by passively 
evolving dust-free elliptical galaxies. Totani et al. (2001b) ar- 
gued that the colors, magnitudes, and counts of these enigmatic 
objects can be best explained by forming elUptical galaxies red- 
dened by screen-Uke dust. It was also found that the surface 
number density of such hyper-extremely red objects (HEROs) 
is roughly the same with that of high-z starbursts discovered by 
SCUBA in sub-mm bands. In the count model of this paper, a 
significant part of the SCUBA sources can be explained by the 
forming elliptical galaxies, and hence our result here quantita- 
tively strengthens the suggestion of Totani et al. (2001b) that 
the SCUBA sources and HEROs have the same origin of form- 
ing primordial elUptical galaxies in the dusty starburst phase. 

6.2. The cosmic star formation history 

Figure 18 shows the cosmic history of global star formation 
rate as a function of redshift, for our model as weU as observa- 
tional estimates coming from Ha, UV, and submilUmeter lumi- 
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nosity density evolution. The elliptical galaxies in our model 
have a strong peak of SFR at their formation redshift, while the 
formation redshift should be distributed in some redshift range 
in reality. Therefore we plot the mean SFR for elliptical galax- 
ies when they are assumed to be formed before z = 2 or 3, as 
indicated by horizontal lines. 

ElUptical galaxies in our model have much higher global SFR 
at z > 3 than the peak at z ~ 1 suggested by the optical observa- 
tions (Lilly et al. 1996; Connolly et al. 1997; Madau, Pozzetti, 
& Dickinson 1998; Steidel et al. 1999; Cowie, Songaila, & 
Barger 1999; Thompson et al. 2001), and the model is still 
consistent with the FIR background detected by COBEMRAS. 
This is in very sharp contrast to a number of previous studies 
(e.g., Gispert, Lagache, & Puget 2000; Takeuchi et al. 2001a; 
Franceschini et al. 2001) which claimed that such high SFR 
beyond z ~ 1 is disfavored from the COBE/FIRAS data. This 
is because the dust temperature of forming elliptical galaxies is 
much higher than expected from the local LiR-Tdust relation, as 
explained in §5.3. The SFR expected from elliptical galaxies is 
close to those inferred from galaxies detected in submillimeter 
observations (Hughes et al. 1998; Barger, Cowie, & Richards 
2000), and this again suggests that a considerable part of form- 
ing elliptical galaxies has been detected by submillimeter ob- 
servations. [Note, however, that the higher points of Barger et 
al. (2000) are corrected by a factor of ^ 1 1 for incomplete- 
ness.] On the other hand, SFR expected from spiral galaxies 
is in reasonable agreement with those obtained by optical (i.e., 
UV in the restframe) observations. Therefore a natural interpre- 
tation of these results is that the cosmic star formation history 
inferred from UV Ught is that for steady or modest star for- 
mation in disks of spiral galaxies, while starbursts in elliptical 
or spheroidal components of galaxies emit their energy mostly 
in the form of dust emission observed in submillimeter wave- 
lengths. 

6.3. On the origin of optical size-luminosity relation 

The key of our success to reproduce the scaling properties 
of local infrared galaxies (i.e., optical/infrared luminosity ra- 
tio, LiR-7dust relation, and infrared luminosity function) is the 
relation between the optical luminosity and size of galaxies, 
which is used as an input in our model. Therefore the read- 
ers may wonder what is the origin of this relation. In fact, this 
scaling relation can be explained by the framework of the stan- 
dard structure formation theory in the CDM universe. One can 
predict the virial size of cosmological objects when the mass 
and formation redshift are specified. Then the observed scal- 
ing relation between Lb and is easily reproduced when some 
reasonable assumptions are made such as a roughly constant 
mass-to-luminosity ratio, proportionality between the virial ra- 
dius and disk effective radius (equivalent to the constant angular 
momentum), and formation of disks at z < 1 (e.g., Mao & Mo 
1998). The scatter along the mean relation can also be repro- 
duced by the dispersion of angular momentum of dark haloes. 
This means that the local scaling relations of infrared galaxies 
also have their origin in the structure formation in the CDM uni- 
verse. This also explains the success of a semi-analytic model 
of Granato et al. (2000) to reproduce the local infrared proper- 
ties. 

6.4. Implications for the TeV gamma-ray observations 

TeV gamma-ray observations provide a unique constraint on 
CBR because very high energy gamma-rays beyond ~TeV are 



absorbed by interaction with the infrared CBR photons to pro- 
duce electron-positron pairs in intergalactic field (e.g., Gould 
& Schreder 1967; Stecker & de Jager 1997, 1998; Salamon 
& Stecker 1998; Primack et al. 1999). Gamma-rays with en- 
ergy e-y mainly interact with infrared photons with wavelength 
A ~ hpe-y/i4m^c^) ~ 1.24e-y^TeV where hp is the Planck 
constant and e^ jev = e7/TeV. Therefore detections of gamma- 
rays at ~ lOTeV for TeV blazars Mrk 421 (z = 0.031) and Mrk 
501 (z = 0.034) (Krennrich, et al. 1999; Aharonian et al. 1999) 
give constraints on the MIR CBR at about 10-20 /im. The thick 
crosses plotted in Fig. 13 are taken from Biller et al. (1998) 
and Renault et al. (2001), and our baseline model is consistent 
with this limit. The CBR flux at 15/im in our baseline model is 
2.5nW/m^/sr, and the lower bound set by the count integration 
in this band accounts for about 80 % of the total CBR of the 
model. 

It should be noted that the model using the empirical LiR-Tjust 
relation predicts higher MIR CBR flux, while the FIR peak is 
lower than the baseline model (Fig. 17). This difference origi- 
nates from the dependence of MIR/FTR flux ratio of the model 
infrared SEDs on the dust temperature. As can be seen in Fig. 1, 
MIR/FIR flux ratio decreases with increasing dust temperature, 
due to the diminishing component of small grain dust and PAH 
features. This trend is physically reasonable as explained in 
§3.2. The difference of our baseline model and that using the 
empirical LiR-Tjust relation is that the dust temperature of star- 
bursting elliptical galaxies in the former is much higher than in 
the latter. Therefore, in the baseline model, the strong thermal 
dust emission from forming eUiptical galaxies has a largest con- 
tribution to the peak of FIR CBR, while it hardly contributes to 
the MIR CBR (Fig. 13). This is why we could reproduce simul- 
taneously the high FIR peak and low MIR CBR flux which are 
fuUy consistent with observations of COBE and TeV gaimna- 
rays. 

On the other hand, the dust temperature does not become as 
high as our baseline model in models employing the empirical 
^iR-7dust relation, and hence MlR/FlR flux ratio of CBR is not 
low enough to match both TeV limits in MIR and tiie COBE 
data in FIR. This prima facie conflict has been noticed in pre- 
vious work based on the empirical LiR-Tjust relations of local 
infrared galaxies (e.g., Stecker 2000), and several possibilities 
have been proposed to solve this conflict, including new physics 
such as a possibility of the broken Lorentz invariance. How- 
ever, our model suggests that the MIR and FIR CBR are made 
by completely different populations of galaxies: low-redshift 
spiral galaxies for MIR while high-redshift starbursting ellipti- 
cal galaxies for FIR, as seen in Fig. 13. Thus in our model the 
conflict is naturally resolved. 

We calculated the optical depth of TeV gamma-rays as a 
function of source redshift (zj) and gamma-ray energy observed 
at z = 0, based on the formulations given in Salamon & Stecker 
(1998). The result is shown in Fig. 19, along with the previ- 
ous calculations of Stecker & de Jager (1998, hereafter SD98) 
based on the empirical LiR-Tdust relation. The baseline model is 
used for the flux and its evolution of the CBR from optical to 
FIR. As expected, the optical depth of our model (indicated as 
TTTTT in the figure from the initials of the authors) is relatively 
smaller at gamma-ray energies of ^ 1-10 TeV where the main 
target photons are in the MIR band, while it becomes larger for 
higher energy gamma-rays interacting with FIR photons, when 
compared with the SD98 calculation. The different two models 
might be discriminated by more accurate measurement of the 
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TeV spectrum of Mrk 421 and/or 501 in future observations. 
Following SD98, we give the result of polynomial fits to the 
optical depth t(zs, e-y), as 

TiZs,ej)= ^ ^ aijilogiQej^TevYQog^oZsy , (14) 

0<K3 0<;<2 

where the values of aij are presented in Table 1 . This fitting for- 
mula gives a reasonable approximation (< 10% at Zs = 0.03 and 

1 < e7,TeV < 20, and < 50% within ranges of 0.03 < Zs < 0.3 
and 0.1 < e-y TeV < 50). More accurate table of the optical depth 
with wider ranges of Zs and is available on request to the au- 
thors. 

7. SUMMARY AND CONCLUSIONS 

In this paper we developed a new model of counts and cos- 
mic background radiation (CBR) of infrared galaxies observed 
by emission from heated interstellar dust, by extending a model 
for optical/NIR galaxies (Totani & Yoshii 2000). Five morpho- 
logical types of galaxies (E/SO, Sab, Sbc, Scd, and Sdm) are 
taken into account and their number densities are normalized 
by type-dependent fi-band luminosity function at the local uni- 
verse. Their luminosity evolution is traced backwards in time 
based on star formation histories inferred from the present-day 
optical/NIR SEDs and chemical properties. Formation epoch 
of galaxies is a parameter for which we tried a redshift range of 

2 <zf < 5. Pure luminosity evolution without number evolu- 
tion is assumed in our baseUne model, but some number evolu- 
tion is also tested in a phenomenological way. The model has 
already been compared comprehensively with the counts and 
redshift distributions of galaxies observed in optical and NIR 
bands, and found to be in reasonable agreement with the data 
(Totani & Yoshii 2000; Totani et al. 2001c). 

Relatively rare populations of AGNs and ULlRGs seen in 
the local universe, whose contribution to counts and CBR is ex- 
pected to be small, are not included in our model, and we tried 
to explain the bulk of infrared data by normal galaxy popula- 
tions at the local universe and their ancestors at high redshifts. 
On the other hand, our analysis strongly indicates that the pri- 
mordial elliptical galaxies are very similar to dusty starburst 
galaxies or ULIRGs with very high dust temperature. 

There are two important new characteristics in our model 
that are different from previous models of infrared galaxies: (1) 
mass scale dependence of dust extinction is introduced by the 
observed size-luminosity relation of optical galaxies, and (2) 
dust temperature is determined by physical consideration of en- 
ergy balance, rather than using the empirical relation between 
the total infrared luminosity (Lir) and characteristic dust tem- 
perature (Tdust) of local galaxies that has been used in a number 
of previous models. As a result, the local properties of infrared 
galaxies, such as optical/infrared luminosity ratios, correlation 
between infrared luminosity and dust temperature, and infrared 
luminosity function, are outputs that should be compared with 
observed data. Indeed we found that our model quantitatively 
reproduces the observed infrared properties at the local uni- 
verse. The key to understand these scahng properties is the 
size-luminosity relation of galaxies; surface brightness and dust 
column density increase with increasing optical galaxy lumi- 
nosity, and hence more massive galaxies should be more dusty. 
This gives a quantitative explanation for the observed correla- 
tion between optical and infrared luminosities. Furthermore, 
massive galaxies should emit more energy as dust emission per 
unit mass of dust, and hence the energy balance inevitably re- 
sults in higher temperature for larger galaxies, as observed for 



local infrared galaxies. The scatter along the mean Ls-rg rela- 
tion is comparable with those in Ltr-L^ and LiR-7dust relations. 
These effects result in much faster increase of 60/xm luminosity 
when Lb is increased, and giving an explanation for the much 
broader shape of 60/im luminosity function than the Schechter 
function of the optical luminosity function. 

Then we predicted faint source counts and CBR composed 
of high-z galaxies. Our baseline model assumes a cosmolog- 
ical model with (/!,f7o,f^A) = (0.7, 0.2, 0.8), pure luminosity 
evolution after formation at zf = 3, and screen distribution of 
dust. We found that this baseline model is in reasonable agree- 
ment with all available data of galaxy counts in six wavebands 
(15, 60, 90, 170, 450, and S50^J.m) and CBR. Therefore our 
model, based only on present-day normal galaxy populations 
and their evolution, reasonably fits to all available data from 
optical to submillimeter wavebands, though some modest num- 
ber evolution may be required for late-type galaxies. The high-z 
starburst galaxies discovered by SCUBA are quantitatively well 
explained in our model by the emergence of starbursts at the 
formation of present-day elliptical galaxies at z ~ 3. 

We also tested that slab-type dust distribution as well as the 
screen dust used in the baseline model, and found that the 
screen model gives much better fit to the observed data of the 
local infrared luminosity function, galaxy counts, and CBR. Al- 
though a pure screen distribution may seem unlikely, it is rather 
reasonable to expect that a part of dust particles behave like 
an effective screen, because of stellar/galactic wind or inhomo- 
geneity of interstellar medium. If it is the case, and the dust 
opacity to optical/NIR is much larger than the unity, the screen 
model should be better than the slab-type prescription. It is 
also supported by extremely red colors of local starburst galax- 
ies (Calzetti, Kinney, & Storchi-Bergmann 1994) or recently 
discovered hyper extremely red objects (Totani et al. 2001b), 
which cannot be explained simply by slab-type dust. There- 
fore we conclude that screen dust is a better phenomenological 
description than slab dust at least in a study of this kind. 

The most drastic difference of our model from previous ones 
is that the dust temperature of starbursting elliptical galaxies 
is predicted to be much higher 40-80K) than that extrapo- 
lated by the empirical LiR-Fdust relation of local infrared galax- 
ies. This is because starbursting elliptical galaxies should emit 
much larger amount of energy as dust emission per unit dust 
mass, than local galaxy populations. On the other hand, such 
high temperature is similar to those found in local ULIRGs or 
high-z dust starbursts observed in sub-mm bands. Thus, our 
result gives a further support to an idea that the progenitor of 
present-day elhptical galaxies or bulges are dusty starbursts. 

There is an important implication for the cosmic star for- 
mation history, which is very different from previous results. 
A number of papers based on the empirical LiR-Tdust relation 
claimed that cosmic star formation rate beyond z > 1 must turn 
over and keep constant or dechne, otherwise it would produce 
too much submillimeter CBR compared with the COBE/FIRAS 
data, by the redshifted dust emission from high-z galaxies 
(e.g., Gispert, Lagache, & Puget 2000; Takeuchi et al. 2001a; 
Franceschini et al. 2001). However, in our model the dust tem- 
perature of forming elliptical galaxies is much higher than in the 
previous models. As a result, although our model assumes very 
strong starbursts at zf ^ 3 in primordial elliptical galaxies and 
the cosmic SFR at zf > 2-3 is even higher than the peak at z ~ 
1 suggested by optical observations, the prediction is in good 
agreement not only with the COBE^RAS CBR measurements 
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but also with the counts and redshift estimation for submillime- 
ter sources revealed by SCUBA. In our model the dusty star- 
bursts in primordial elliptical galaxies hardly contribute to the 
cosmic SFR measured by optical observations. 

Another result that is significantly different from previous 
studies is the smaller ratio of MlR/FlR CBR flux compared 
with models based on the empirical Z-m-Tdust relation of local 
infrared galaxies (e.g., Stecker 2000). It should be noted that 
our model does include the warm/nonequilibrium components 
of dust emission that dominate the thermal emission from big 
grain dust at the MIR region, while several previous models 
gave very low MIR CBR flux as a valley in the CBR spec- 
trum simply because they did not take them into account (e.g., 
MacMinn & Primack 1996; FaU, Chariot, & Pei 1996). The rea- 
son we get this result is the trend that the warm/nonequilibrium 
components become less significant compared with thermal 



big-grain dust emission, with increasing temperature of dust. 
This is inferred from the infrared spectrum of local galaxies, 
and it is also expected from physical consideration. If this 
trend is taken for granted for high-z galaxies, the MIR/FIR flux 
ratio of dust emission from forming elliptical galaxies should 
be much smaller than that of local infrared galaxies, since the 
dust temperature of starbursting eUiptical galaxies is found to 
be very high in our model. Therefore these galaxies dominate 
in the FIR peak of CBR while they have very small contribu- 
tion in MIR CBR. This effect is significant enough to resolve 
the prima facie conflict (Stecker 2000) between the upper hmits 
on MIR CBR from TeV gamma-ray observations and FIR CBR 
detections by COBE/DIRBE and FIRAS. 

The authors have financially been supported by the JSPS Fel- 
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TABLE 1 
POLYNOMIAL FITS TO T{ej,Zs) 



J 




aij 


a2j 


fl3j 





1.11 


0.85 


-0.28 


0.32 


1 


1.16 


0.42 


0.87 


-0.53 


2 


0.04 


0.21 


0.35 


-0.23 



Note. — This fitting formula is accurate within 10% at Zs = 0.03 and 1 < e-y < 20, and within ~ 50% for ranges of 0.03 < Zs < 0.3 
and 0.1 <e.y<50TeV. 
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Fig. 1 . — The infrared .spectral energy distribution (SED) of model galaxies, 
for various values of the temperature of big grain dust in thermal equilibrium 
(indicated in the figure). The nonthermal emission coming from small grain 
dust and PAH features are also taken into account. 
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Fig. 2. — Infrared spectral energy distribution (SED) of the Milky Way (top 
panel) and M82 (bottom panel). The model curves are fit by our model SED 
described by a characteristic temperature of cold, big grain dust (indicated in 
the panels), while the nonthermal emission coming from small grain dust and 
PAH features are also added. The observed data are from Dwek et al. (1997) 
for the Milky Way, and Thuma et al. (2000) for M82. 
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Fig. 3. — The same as Fig. 2, but for other four nearby infrared galaxies. 
The observed data are taken from Klaas et al. (1997) for Arp 220 and Arp 244, 
Prayer et al. (1999) for VV 114, and Yun et al. (2000) for BR 1202-0725, 
respectively. 
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Fig. 4. — Correlation plot of the bolometric infrared luminosity versus opti- 
cal luminosity (B-band) of galaxies. The data are from Saunders et al. (2000). 
The dot-dashed line is that corresponds to Li^/vL^{B) — 1. The thick solid 
and dashed lines are the prediction of our model assuming the screen and slab 
distribution of dust, respectively. The thin solid and dashed lines are the same 
with the thick Unes, but in the cases where the size of galaxies are shifted by 
itliT of observed dispersion in logr^. 




1 I I I I I 1 ' I 1 t I 1 1 I 1 r 1 1 I t 1 1 I I 

108 109 IQIO 10" 1012 

Bolometric L|p [h"^ Lg] 

Fig. 5. — Correlation plot of the bolometric infrared Itmrinosity versus tem- 
perature of cold, big grain dust estimated from 60 and 100 /im fluxes. The 
data are taken from Saunders et al. (2000). The thick solid and dashed lines 
are the prediction of our model assuming the screen and slab distribution of 
dust, respectively. The thin solid and dashed lines are the same with the thick 
Unes, but in the cases where the size of galaxies are shifted by ±1it of ob- 
served dispersion in logr^. The dot-dashed line is a simple linear fit to the 
data, which is used to determine the dust temperature based on the empirical 
i-iR-Tdust relation (eq. 13). 
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Fig. 6. — Infrared luminosity function of galaxies at 60/im. The data are 
taken from Saunders et al. (1990). The solid and dashed lines are the pre- 
diction of our model assiuning the screen and slab distribution of dust, re- 
spectively. Considering the uncertainty in the normalization of luminosity 
function, the normalization of model curves are set to match the data in the 
range of i/L^(60/im) < 10"ft~^LQ. [More Ivmiinous galaxies are removed 
from normaUzation fit to avoid contamination of AGNs which is expected to 
be significant for ultraluminous galaxies (Veilleux, Kim, & Sanders 1999).] 
The dot-dashed line is the prediction from the Schechter-type B-band luminos- 
ity function when the infrared luminosity is simply proportional to the B-band 
luminosity in all galaxies as vLtj(B)/vLv((f)fxm)=l. 
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Fig. 7. — The evolution of infrared luminosity and characteristic temperature 
of big grain dust predicted by our model as a function of time after the forma- 
tion, for five types of galaxies: E/SO (solid), Sab (short-dashed), Sbc (long- 
dashed), Scd (dot-short-dashed), and Sdm (dot-long-dashed). The upper-left 
panel shows the evolution of total infrared luminosity of our model assuming 
the screen-type dust. The upper-right panel shows the temperature evolution 
predicted when one uses the empirical Lir-Tjusi relation observed for local in- 
frared galaxies. The lower-left and -right panels show the predictions of dust 
temperature evolution based on our model taking into account the energy bal- 
ance between absorption and re-emission of stellar light, for the screen- and 
slab-type dust distributions, respectively. All models are for galaxies whose 
luminosity is Mb = -20 at z = 0. 
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Fig. 8. — Faint galaxy counts in infrared and submillimeter bands. Tlie solid line is the prediction by our baseline model with a cosmological model of 
(A, Qqi^a) = (0.7,0.2.0.8), pure luminosity evolution, formation redshift of Zf = 3, and screen-type distribution of dust. The solid line is the total of all galaxy 
types, while the other five lines are for individual types: E/SO (dotted). Sab (short-dashed), Sbc (dot-dashed), Scd (three-dot-dashed) and Sdm (long-dashed). The 
data are taken from the following references: Rush et al. (1993), Oliver et al. (1997), Flores et al. (1999a, b), Clements et al. (1999), Aussel et al. (1999), Altieri 
et al. (1999), Elbaz et al. (1999), and Serjeant et al. (2000) (15 ^tm); Rowan-Robinson et al. (1991) (60 fim); Juvela et al. (2000), Matsuhara et al. (2000), Efstathiou 
et al. (2000), and IRAS 100 ^m galaxy counts (90 fim); Stickel et al. (1998), Kawara et al. (1998), Puget et al. (1999), Juvela et al. (2000), Matsuhara et al. (2000), 
and Dole et al. (2001) (170 ^m); Small et al. (1997), Barger, Cowie, & Sanders (1999), and Blain et al. (2000) (450 ^m); and Smail, Ivison, & Blain (1997), Hughes 
et al. (1998), Blain et al. (1999), Bales et al. (1999, 2000), Holland et al. (1999), and Barger et al. (1998, 1999) (850 fim). 
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Fig. 9. — The same as Fig. 8, but for different cosmological models: the A-dominated flat universe with {h,Uo,^A) = (0.7, 0.2, 0.8) (sohd, the baseline model 
shown in Fig. 8), an open universe with (0.7, 0.2, 0.0) (dotted), and an Einstein-de Sitter universe with (0.5, 1, 0) (dashed). 
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Fig. 10. — The same as Fig. 8, but for different formation redshift of galaxies: zf = 2 (dotted), zf = 3 (soUd, the baseline model shown in Fig. 8), and zf = 5 
(dashed). 
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Fig. 1 1 . — The same as Fig. 8, but for the cases of number evolution of galaxies: ?j = (solid, the baseline model without number evolution shown in Fig. 8), 
T) = 1 (dotted), and ri = 2 (dashed). See text for the modeling of number evolution. 



A BRIDGE FROM OPTICAL TO INFRARED GALAXIES 



21 



10 



10 



Co 

A 



lo'- 



10" - 



o lO-' 



g 



5 10"=' - 



' 1 ' 


1 ' 1 ' 

15 /^m 










— std 




emp 


v 

1 ^kv. 


— slab 

- - d-Z 

1 



10 



10 r 



10 



-6 



Flux density 5"^ [Jy] 



10 



10 



Co 
A 



10^ 



° 10' - 



S 10' 



10^ 



' 1 ' 1 


' 1 ' 

90 /^m 






















1 , 1 





10 



10 



10^ 



° 10'- 



5 10' 



60 /um 




Flux density [Jy] 



10^ 



' 1 ' 1 


1 ' 

170 ;tim 














■ v\ 


























\ _ 


1 , 1 


\ 



Flux density [Jy] 



Flux density [Jy] 




Fig. 12. — The same as Fig. 8, but with different model prescriptions from the basehne model (sohd line, and shown in Fig.8): empirical LiR-Tdust relation 
(dotted), slab-type dust (dashed), and non-proportional dust/metal ratio (dot-dashed). 
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Fig. 13. — Spectrum of the cosmic background radiation (CBR) from op- 
tical to FIR bands. The solid line is the prediction by our baseline model, 
and the other five lines show the contribution from individual types of galax- 
ies, as indicated in the figure. References for the observational data are as 
follows; UV data shown by open circles: Martin, Hurwitz, & Bowyer (1991; 
1600 A) and Armand et al. (1994; 2000 A). Ranges estimated from integration 
of optical/NlR galaxy counts: Totani et al. (2()01a, thick error bars without 
symbols). Reports for diffuse CBR detections at optical and NIR bands are 
shown by open diamonds: Bernstein, Freedman, & Madore (2002; 0.3, 0.55, 
and 0.8/im), Gorjian et al. (2000; 2.2 and 3.5/im), Wright & Reese (2000; 2.2 
and 3.5//m), Cambrgsy et al. (2001; 1.25 and 2.2//m), and Wright (2001; 1.25 
and 2.2/im). Lower limits from infrared/submillimeter galaxy count integra- 
tions shown by open squares: Gispert et al. (2000; 6.5 /im), Clements et al. 
(1999; 12 ^tm), Elbaz et al. (1999; 15 ^tm), Lonsdale et al. (1990; 60 /.tm), 
Dwek et al. (1998; 90 ^m), Puget et al. (1999; 175 ^m), Barger et al. (1999; 
850 ^lm). MIR upper limits from TeV gamma-ray observations shown by thick 
crosses: Biller et al. (1998) and Renault et al. (2001). COBE/DIRBE detec- 
tions at 140 and 240 /im: Hauser et al. (1998, filled squares) and Lagache et 
al. (2000, filled triangles). COBB/FIRAS detection at 100 /im - 1 mm shown 
by the shaded area: Fixsen et al. (1998). 
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Fig. 14. — The same as Fig. 13, but for different cosmological models: the 
A-dominated flat universe with (A, Slo, ^^a) = (0.7, 0.2, 0.8) (soUd, the baseline 
model shown in Fig. 13), an open universe with (0.7, 0.2, 0.0) (dotted), and an 
Einstein-de Sitter universe with (0.5, 1, 0) (dashed). 
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Fig. 15. — The same as Fig. 13, but for different formation redsliift of 
galaxies: zf = 2 (dotted), zf = 3 (soUd, the baseUne model shown in Fig. 13), 
and ZF = 5 (dashed). 
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Fig. 17. — The same as Fig. 13, but with different model prescriptions from 
the baseline model (soUd Une, and shown in Fig. 13): empirical ins-rd^st re- 
lation (dotted), slab-type dust (dashed), and non-proportional dust/metal ratio 
(dot-dashed). 
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Fig. 16. — The same as Fig. 13, but for the cases of number evolution of 
galaxies: r; = (solid, the baseline model without number evolution shown 
in Fig. 13), ?? = 1 (dotted), and 77 = 2 (dashed). See text for the modeling of 
number evolution. 
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Fig. 18. — Cosmic history of global star formation rate. The solid curves 
are the prediction of our baseline model for spiral galaxies (i.e., except for el- 
liptical galaxies), with several values of zf- The two horizontal solid lines are 
mean SH? of elliptical galaxies when they are assumed to be formed before 
Zf= 2 or 3, respectively. The observed data for dust-uncorrected SFR from 
Hq or UV luminosity density (i.e., unabsorbed SFR) are shown by open sym- 
bols: Gallego et al. (1995, square), Tresse & Maddox (1998, triangle), Treyer 
et al. (1998. four-edged star), Lilly et al. (1996, open circles), Connolly et 
al. (1997, upside-down triangles), Madau, Pozzetti, & Dickinson (1998, di- 
amonds). Steidel et al. (1999, five-edged stars), Cowie, Songaila, & Barger 
(1999, pentagons). The SFR density inferred from submillimeter observations 
(i.e., hidden SFR) are shown by filled symbols: Hughes et al. (1998, circle); 
Barger, et al. (2000) before (triangles) and after (squares) the completeness 
corrections. The local hidden SFR density from ULlRGs is shown by the 
filled diamond (Barger et al. 2000). These observed data have been corrected 
for the cosmological parameters chosen here [(/i, Hoi^^a) = (0.7,0.2,0.8)]. 
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Fig. 19. — The intergalactic optical depth of very high energy gamma- 
rays to the absorption by interaction with optical/infrared cosmic background 
radiation, as a function of the source redshift (zs) and the gamma-ray energy 
observed at z = 0. The solid lines are the calculation based on our baseUne 
model, with the source redshifts indicated in the figure (indicated as TTTTT 
from the initials of the authors). For comparison, calculations with 'higher 
ERF SED' (dashed line) and 'lower IIRF SED' (dot-dashed line) by Stecker 
& de Jager (1998) are plotted for Zs = 0.03, 0.1, and 0.3 (from the bottom to 
the top). 



